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Abstract. Intense electromagnetic pulses interacting with a plasma can create a 
wake of plasma oscillations. Electrons trapped in such oscillations can be accelerated 
under certain conditions to very high energies. We study the conditions for the 
Wakefield acceleration to produce MeV electrons in atmospheric plasmas. This 
mechanism may explain the origin of MeV or runaway electrons needed in the current 
theories for the production of Terrestrial Gamma ray Flashes. 
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Electromagnetic pulses interacting with a plasma can create a wake of plasma oscillations 
through the action of the nonlinear ponderomotive force [1]. Electrons trapped in 
the wake can be accelerated to high energies. In the laboratory very high-power 
electromagnetic radiation from lasers is used to accelerate electrons to high energies 
in a short distance. However, in atmospheric plasmas the distance range goes from 
metres to kilometres, and the conditions for the wakeheld acceleration takes place are 
different as we will make clear in this work. 

The acceleration of the electrons in atmospheres has been the subject of attention 
since the discovered of Terrestrial Gamma Flashes (TGF’s). TGF’s were hrst discovered 
by accident in 1994 [2] and only limited experimental data is available. One of the goals 
of the European Space Agency (ESA) mission ASIM [3] will be to provide experimental 
data on TGF’s. Satellite measurements suggest that TGF’s typically last from a few 
tens of microseconds to a few milliseconds (50 /rs - 1 ms) [21IH El E] ; exhibit an energy 
spectra extending from 1 MeV up to 100 MeV I211Z1I81E1 [To] , over a half conic-angle 
of 30° — 40° im [12] and have a fluence of 0.1-0.7 photons/cm^. They seems to appear 
around 15 km near the upper regions of thunderclouds uni [HE]. Any theory about 
the TGF generation should be able to successfully explain the observed duration, energy 
spectra, and photon fluence. There is some agreement that the emission is produced 
via Bremsstrahlung HSl H when high energy electrons collide with nuclei in the air 
releasing energy. The satisfactory theory for TGF’s needs to explain the origin of the 
high-energy electrons [E] which can trigger the ignition of the gamma flashes. One 
explanation is given by the feed-back theory [7| in which seed electrons are provided by 
cosmic rays. Other possibilities are leader/streamer theory [19], in which seed electrons 
are produced in the streamer head and accelerated in the stepped leader electric held, 
and the current pulse version ESI- Here we present an alternative scenario where 
those energetic electrons come from the interaction of a electromagnetic pulse with 
an existing plasma previously created. Although the idea of TGF’s production by an 
electromagnetic pulse due to a lightning return stroke was explored previously IZH, the 
novelty of our proposal is the interaction of the pulse with a plasma already present in 
the atmosphere and created for example by another discharge or any electromagnetic 
activity. Since the interaccion of the electromagnetic pulse is with a plasma, under 
certains conditions the braking in the motion of the electrons by friction does not play 
a role, as we will show in this work. 

The outline of the paper is the following. First we revise how an electromagnetic 
wave can create a wakeheld in a plasma. We study the conditions under which the 
acceleration of electrons to MeV energies is possible in the presence of atmospheric 
plasmas. The acceleration mechanism is demonstrated through computer simulations. 
We present the scenario for the production of the plasma and the electromagnetic pulse 
in planetary atmospheres, and in particular for the Earth atmosphere we give some 
predictions. Finally we summarise the main results and end with some conclusions. 


Wakefield acceleration in atmospheric plasmas 

2. Conditions for electromagnetic propagation of waves into a plasma 


3 


Here we study the interaction of an electromagnetic wave with a plasma. Let us consider 
how a plasma modihes the propagation of electromagnetic waves [22]. The starting point 
to describe the evolution of a plasma is the Vlasov equation, which up to hrst order in 
the expansion parameter 1/Nd, reads 
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being /j(r, v,t) the distribution function of the j species. The right hand side of ([^ 
represents the collisional terms. The parameter Nd is the number of particles in the 
Debye sphere. For electrons, Nd = 47rneA|,g/3, where Ue is the electron density and Xoe 
the electron Debye length. An useful expression of the Debye length under equilibrium 
conditions is given by A^e = bQ-^T/ng (m), where T is the equilibrium temperature of 
the electrons in Kelvin and Ug is expressed in m“^ [23], so we have the condition 


Nd ~ 1.38 X 10^ > 1, 


( 2 ) 


for Q to be valid. We recall that the condition to have a plasma is Xdc ^ L, being L 
the characteristic dimension of the system. 

From Q it is a standard procedure to calculate the hrst moment equations. We 
will assume that collisions do not change the number of species, so in the averaging we 
will take 
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Hence, the change of the momentum of the j species becomes 
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where nj is the density and Uj the velocity of the j species and pj the pressure. 

Let us consider a plasma composed of ions with positive charged +eZ and electrons. 
We will assumed that ions form a hxed background with density n, = n^jZ, where 
is the corresponding electron density. Then we only need to treat the dynamics of the 
electron huid. We can investigate the damping of an electromagnetic wave of the form 
E(r) exp{—iut) considering the linearised plasma response. Writing 


d 

Wf 

where Vgi is the collisional frequency of the scattering of electrons by ions, from 
first order we get 
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The solution of Q results 
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The plasma conductivity a can be calculated from the current density of the plasma 
j = -erieUe- Using 0 


oo. 


J = ^£o 


pe 


UJ + iUe 


E(r)e-*"'* = aE(r,f), 
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where the plasma frequency is defined as 
ule = neeV(£0"ie), 


(7) 


and the plasma conductivity cr is a complex quantity, 

a = i£Q —^—. (8) 

UJ + lh>ei 

From Maxwell’s equations for harmonic fields, 


V X E(r) = -ia;B(r), 


V X B(r) = /iocrE — i-^E(r). 

c2 

(9) 

The second equation in 0 can be written as 


V X B(r) = -i^£E(r), 

(10) 

where 


£ = l ^ 

UJ{UJ + Wei) ’ 

(11) 

is the dielectric function of the plasma. Taking the curl 

in equations 0 and (10), one 


gets the wave equations 
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V"E(r) - V (V ■ E(r)) + ^£E(r) = 0, 
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V^B(r) + -Ve X (V X B(r)) H —^£B(r) = 0, 


( 12 ) 


that give the spatial behaviour of the electric and magnetic fields in the plasma. 

In the case of a neutral plasma with an uniform density Ve = 0 and V ■ E = 0, for 
harmonic electromagnetic waves E(r) ~ exp(fk-r), the equations (12) yield uj'^e = 

So using 

= (13) 

V UJ / 

in which it is assumed that Vei ^ Equation (13) means that electromagnetic waves 
are damped. Writing uj = uj^ — iv 12^ where v is the damping rate, equation (13) becomes 

ujr = + A:2c2, 

(14) 
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The damping rate can be computed from the zero-order distribution function. This 
allows to get an expression for the collision frequency, namely [221 [23] 
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where we use Ve = ^/KsT/me for the thermal velocity of electrons. The factor A is 
the ratio of the maximum and minimum impact parameters. The maximum impact 
parameter rmax is given by the Debye length \oe, as the Coulomb potential is shielded 
out over that distance. The minimum impact parameter is given by the classical distance 
of closest approach = Ze’^/A'neomf.vl, which averaging over all the particle velocities 
and assuming a Maxwellian distribution yields f^m = Ze^/ (WtieoKeT), so 


A = 


A 


De 




(16) 


For typical values of plasmas, InA 10. 

We must notice that (15) is only an approximation where the zero order electron 
distribution is taken Maxwellian. For non equilibrium processes, the collisional damping 
rate could be less. For example in the case of a super-Gaussian distribution, the 
collisional damping is reduced by a factor of 2. 


3. The acceleration of electrons in atmospheric plasmas 


In this section we will study how electrons in a plasma are accelerated by an 
electromagnetic pulse propagating through the plasma. We will first set the equations 
to describe the electron dynamics and then compare their predictions with some particle 
in cell code simulations to verify the results. 

Let us consider an electromagnetic pulse which is propagating in the x-direction 
inside a plasma, 

E(r, t)= [ dkA{k)e^^^^-^^^ey + c.c., (17) 


where A{—k) = A{k)* and A{k) for A; > 0 is nonzero only in the vicinity of a 
central wavenumber k^. From the relation (14) we take for the pulse a lead frequency 
ojq = (Wpg + Expanding around ko yields the following expressions for the 

electromagnetic field. 


E(r, t) = Eq f{x - Vgt) cos[(i:no - Vgko)t] By, 

B(r, t) = f{x - Vgt) cos[(a;o - Vgko)t] e^, (18) 


where 


Vg = I 


(-) 


C^kn 
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(19) 


k=kg 

is the propagation speed inside the plasma, £ is the plasma dielectric function ( [II| , and 
j{x) = Eyij, 0)/Eo gives the shape of the pulse, Eq being the amplitude of the electric 
field. 


As discussed in the previous section, the condition ujq > Upe is required for the 
propagation of the electromagnetic pulse inside the plasma. In the following we will 
consider ujq > 10 ujpe. For loq = lOojpe, Vg differs from the speed of light in vacuum 
in less than 1%. For larger values of (Uq, the pulse propagation speed is even closer to 
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c. The electromagnetic pulse will start interacting with the electrons at t = 0. The 
interaction will hnish at a time of the order of 2n/uo. Since Uq — Vgko = the 

cosine term in (18), cos[(u;o — Vgko)t] = cos[(ci;pe/i^o)^i^o^] ~ 1 can be safely neglected 
throughout the process. We will also assume that the damping rate (15) of the pulse is 
negligible. Those assumptions will impose some constrains that we will discuss in the 
next section. With these approximations, the electromagnetic held can be written as 

E{r, t) = Eof{x - ct)ey, (20) 

ct)e^. (21) 


B(r,t) = —f{x 
c 


Since the pulse frequency is larger than the plasma frequency, the dynamics of the 
electrons will be collisionless while interacting with the pulse, so can be modelled by 


dt 


ev X B, (22) 

where p = rueyv is the relativistic momentum of the electron and 7 = (1 — 

Equation (22) can be written as 


dv ^ (I ^ 

dt rUe 
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We dehne dimensionless variables t = teEf)/{m^c) and f = ycEq/{ me(?)-, so that the 
scaled velocity is just v = v/c. Using (20) and (21), the equation of motion (23) 
becomes 


dt 
dv. 


v'^ivn, 


VxVy 




= -Vl -u2(l - Uy - Vx) f[u:{x - t)], 

= -x/1 

dt 

where f{kox) = f{x) takes into account the assumed form of the wave-packet, and 

meCOjQ 


bJ = 


eEn 


(24) 


(26) 


We are interested in pulses which are able to accelerate electrons to kinetic energies 
K = meC^[(l — — 1] in the order of several MeV after the interaction time, which 

is T ^ 2njCj. Therefore, in this model the interaction with the pulse is controlled by 
the dimensionless parameter 0. Pulses with different amplitudes and frequencies but 
same ratio uq/Eq (and shape) produce the same acceleration. 

Figure [ 1 ] shows the relativistic kinetic energy obtained for a free electron initially 
at rest at the origin accelerated by a pulse with the shape of the form 


f{ux) = cos[c5(a; - f/A)][H{x + T/2) - H{x)], 
which is unipolar (see right of Fig. [^, and also for the bipolar shape 
f{ux) = cos[(h(x - f/4:)][H{x + f) - H{x)], 


(26) 


(27) 
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Figure 1. Kinetic energy vs time of an electron accelerated by a unipolar (dotted 
lines) and a bipolar (crosses) EM pulse, for several oj. The pulses, depicted in the 
right side, are assumed to be translation invariant in the directions perpendicular to 
the direction of propagation (along the x-axis). 

X (m) 


36.72 36.75 36.78 36.81 



Figure 2. Energy of electrons in a plasma of density no = 10^° m“^ after an 
acceleration time interval of 1.44 x 10“^ s due to an electromagnetic pulse wq = 
27 r X 10^ rad/s. The continuous line corresponds to a unipolar pulse of amplitude 
Eq = 10®V/m and its abscissa is the upper one. The dashed line correspond 
to Eq = 10® V/m, being its abscissa the lower one. The plasma frequency is 
ujpe = 5.6 X 10® rad/s. 


where H{x) is the Heaviside step function. The results were obtained by solving 
numerically (24). The simulations show that for values O < 0.1, regardless whether 
the shape is unipolar or bipolar, the kinetic energy keeps growing after five periods Tq, 
which is an indication that the electron is actually trapped by the pulse, soon achieving 
very high energies. For a; ~ 1, the electromagnetic pulse is not so strong as to carry 
away the electron. It is still able to produce accelerations to energies in the MeV 
scale, though for the bipolar pulse the electron acceleration is later reversed by the own 
electromagnetic packet. For larger values of the reduced frequency, a; > 10, the impulse 


produced by the pulse is well below the MeV range. 

These results were numerically verihed using the Finite Difference Time Domain 
Particle in Cell software VORPAL |2S]. Figure [^shows the kinetic energy prohle after a 
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time interval of 1.44 x 10“^ s, for a plasma of density no = 10^° m“^ with a nnipolar pnlse 
defined by the bonndary condition E(0, y, z, t) = Eq sin(a;ot) H{n/uJo — t) By, where H is 
the Heaviside step fnnction. For Eq = 10^ V/m and O = 1.071, as predicted, the pulse is 
able to produce electrons just in the MeV scale, whereas the pulse with Eo = 10® V/m 
yields a) = 0.1071, able to accelerate free electrons to very large energies. Note that for 
the case of the higher electric field amplitude, the acceleration occurs at larger distances. 

4. A proposed scenario for Gamma flashes in planetary atmospheres 

We have shown in previous section how an electromagnetic pulse, interacting with a 
created plasma, is able to accelerate electrons to high MeV energies provided some 
conditions are fulhlled. We propose the following scenario to this process takes place in 
planetary atmospheres. 

Some electromagnetic activity has to produce the pulse and the plasma, for example 
lightning discharges cloud to cloud, or cloud to earth. When the pulse is created, it 
propagates at basically the speed of light through the planetary atmosphere. If that 
pulse hnds in its way a plasma, then the wakefield mechanism might be fired. We 
can speculate that the origin of the plasma is another discharge, or even the same 
discharge that creates the pulse. We will come back to this point in the next section 
when studying the Earth atmosphere case. The high energy electrons accelerated by 
the wakefield, leaving the plasma region, can interact with neutral molecules, creating 
an additional shower of secondary electrons, other particles and gamma radiation via 
Bremsstrahlung. In Figure a cartoon picture summarises the whole cascade of events, 
leading eventually to the creation of a gamma flash. 

From previous sections we need the following conditions for this scenario to work: 

• The numbers of particles in the Debye sphere must be large, 

Nd > 1. (28) 

• The Debye length must be much smaller than the characteristic size of the plasma 
region, 

^De L. (29) 

• The frequency of electron and ions collisions must be smaller than the 
electromagnetic pulse angular frequency, 

Wi < ^^ 0 - (30) 

• The damping rate must be negligible, 

/ >2 

n = -^Uie < 1 . 

^0 

• For the acceleration of electrons to MeV energies, we need the condition 
meCOJo 


( 31 ) 
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Figure 3. Cartoon outlining the events leading to the creation of a TGF, courtesy of 
D. Burton. 


5. Predictions for TGF’s 


In this section we consider the ignition of TGF’s in the Earth atmosphere by the 
wakeheld mechanism. The occurrence of transient luminous events (TLE’s) are 
candidates for the electromagnetic pulse source. At the troposphere altitudes, from 
0 to 20 km, we have conventional cloud to ground lightning discharges of both types, 
positive and negative ones, and cloud to cloud discharges. Closer to the stratosphere, 
from 20 to 40 km, there are blue jets [26]. At higher altitudes, at the mesosphere, we 
hud sprites and, at the ionosphere, elves [27] . 

Due to the fact that it seems that TGF’s appear at 15 km as pointed at the 
introduction, we will assume that the electromagnetic pulse is produced by lightning 
activity. We will take for the peak intensity of the pulse 5 x 10® V/m [28]. From 
relation (32), we can get the condition wq < 3 x 10® rad/s. According to ^9] |30], the 
range of average frequency spectra of return strokes is 0.2 — 20 MHz, so we can have a 
pulse with fo = ujo/2n < 5 X 10^ Hz. 

To fulfill (30) and (31), as in the simulations, we take cOpe = ujq/W. It turns out 
that there is a critical maximum electron density of the plasma Uec ~ 3 x 10^^ m“® which 
allows wakeheld acceleration. With those values, InA = 14, where we have taken the 
mean atomic number of the air Z = 7.2, and the electrom plasma temperature 10^ K 


The characteristic plasma length is = 1 cm and (30) = 100 rad/s which is 

much less than cjq so collisions can be negleted. 

Gan a plasma of of size bigger than 1 cm and electron density less than 3 x 10^^ m“® 
be created? At altitudes around 15 km, the air density in the Earth is bigger than 
3 X 10^^ m“®, so the total ionisation of the air will give a too dense plasma. However, 
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the rapid expansion of a plasma, due to the shock wave created by a lightning discharge 
could account for such a low density plasma. At higher altitudes, the plasma density 
can be lower, but the electromagnetic pulse normally decays inversely with the distance 
as 1/r. Then if the electromagnetic pulse is created at large distances from the plasma 
region, the held will not be able to trigger the wakeheld acceleration when reaching the 
plasma. We have considered only the helds created by cloud to ground dischages in 
the discussion, but there are other sources of TLE’s which might provided the necessary 
pulses. At the present state of experimental data, we cannot go further in this discussion. 

6. Further discussions and conclusions 

In this paper we have studied the conditions for the creation of MeV electrons in 
atmospherics plasmas. An intense electromagnetic pulse interacting with the plasma can 
create a wake on the plasma. Electrons trapped in such oscillations can be accelerated 
under certain conditions to high energies. We have shown that those electrons could 
reach energies in the MeV range, thus being able to ignite gamma bursts. We propose 
a scenario for the wakeheld acceleration to take place in planetary atmospheres, where 
the pulse and the plasma are generated by electromagnetic activity. 

Applying our theory to the origin of MeV or runaway electrons for the production of 
Terrestrial Gamma Flashes (TGF’s) we predict that an electromagnetic pulse of 5 x 10^ 
V/m peak can account for the ignition of TGF’s provided that the pulse encounters a 
plasma of low electron density. 

The idea of TGF’s production by an electromagnetic pulse from a lightning return 
stroke propagating into the air was explored previously 1211 . The conclusion was that 
very high electromagnetic held amplitudes are necessary to overcome the friction due 
to collisions with the nonionized gas. 

The novelty of our proposal is the interaction of the pulse with a plasma 
already present in the atmosphere and created for example by another discharge or 
any electromagnetic activity. In order to accelerate electrons to MeV energies, the 
amplitudes of the pulses depend on the density of the plasma. Under the conditions 
studied in this work, the wakeheld mechanism trigger the ignition of TGF’s. We stress 
that if those conditions are fulhlled, the pulse can propagate inside the plasma without 
damping and electrons will be accelerated dragged by the pulse. 

In order to address and clarify the situation, there is scheduled the Atmosphere- 
Space Interactions Monitor (ASIM ESR) mission on the ISS external facilities of the 
Golumbus module. It will be an Earth observation facility devoted to the study of 
severe thunderstorm and their role in Earth’s climate. Among other activities, ASIM 
will study giant electrical discharges (lightning) in the high-altitude atmosphere above 
thunderstorms, and the precise location of TGF’s [3j. 
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